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Hayashi and Carthew (Nature 431 [2004], 647) have shown that the 
packing of cone cells in the Drosophila retina resembles soap bub- 
ble packing, and that changing E- and N-cadherin expression can 
change this packing, as well as cell shape. 

The analogy with bubbles suggests that cell packing is driven by 
surface minimization. We find that this assumption is insufficient to 
model the experimentally observed shapes and packing of the cells 
based on their cadherin expression. We then consider a model in 
which adhesion leads to a surface increase, balanced by cell cortex 
contraction. Using the experimentally observed distributions of E- 
and N-cadherin, we simulate the packing and cell shapes in the wild- 
type eye. Furthermore, by changing only the corresponding param- 
eters, this model can describe the mutants with different numbers 
of cells, or changes in cadherin expression. 

Drosophila retina development | cell shape | surface mechanics 

Cell adhesion molecules are necessary to form a coher- 
ent multicellular organism. They not only hold cells to- 
gether, but differential expression of different types of these 
molecules plays a central role during development. Members 
of the cadherin family are the most widespread molecules that 
mediate adhesion between animal cells, and their role has been 
demonstrated in cell sorting, migration, tumor invasibility, cell 
intercalation, packing of epithelial cells, axon outgrowth and 
many more [U El [3l HI E] [6]. We here focus on the role of 
adhesion in the determination of epithelial cell shape 7 . 

In the compound eye of Drosophila, the basic unit, the om- 
matidium, is repeated approximately 800 times. All omma- 
tidia have the same cell packing, which is essential for correct 
vision. The ommatidium consists of four cone cells, which 
are surrounded by two larger primary pigment cells. These 
'units' are embedded in a hexagonal matrix, constituted by 
secondary and tertiary pigment cells, and bristles (c.f. SI Fig. 

0)- 

Two of us [9] showed that cadherin expression influences 
ommatidial cone cell packing. Two cadherin types, E- and 
N-cadherin, are expressed in different cells: all interfaces bear 
E-cadherin, while N-cadherin is present only at interfaces be- 
tween the four cone cells (SI Fig. [7f. Cadherin-containing 
adherens junctions form a zone close to the apical cell surface, 
allowing the retina epithelium to be treated as a 2D tissue. 
In the wildtype and in i?oi-mutant ommatidia with two to 
six cone cells, these cone cells assume a packing (or topology, 
that is, relative positions of cells) strikingly similar to that of 
a soap bubble cluster. When cadherin expression is changed 
in a few or all of the cells, the topology can change. More 
frequently, only the geometry (individual cell shapes, contact 



angles at the vertices, interface lengths) changes. 

The soap films between bubbles are always under a posi- 
tive tension, 7 > 0. This surface tension describes the energy 
cost of a unit of interface between bubbles, and drives their 
packing. At equilibrium, in a 2D foam layer, soap bubbles 
meet by three at each vertex, since four-bubble vertices are 
unstable [101111] . In addition, since 7 is constant and the same 
for all interfaces, bubble walls meet at equal (i.e. 120°) angles. 
More precisely, the surface energy (or rather the perimeter en- 
ergy, for a 2D foam) is 7P, where P is the total perimeter of 
soap films. The foam reaches equilibrium when it minimises 
P (since 7 is constant), balanced by another constraint fixing 
each bubble's area. 

It has been proposed that cells minimize their surface, 
like soap bubbles [121 1131 114] . Since the surface mechanics 
of bubbles are quite simple, they can easily be described in 
a model. However, calculating the equilibrium shape of a 
cluster of more than four bubbles is difficult [13] : for this pur- 
pose, we use a numerical method [151 116] , in order to test if 
cell patterning is based on surface minimization. Here, the 
only biological ingredient is differential adhesion [l] 1131 117] : 
an interface between two cells has a constant tension, that is 
lower when the adhesion is stronger [151 116] . 

Cells, however, differ greatly from bubbles, both in their 
membrane and internal composition. Surface tension has been 
shown to be determined up to a large extent by the cortical 
cytoskeleton [181 1191 1201 121] . Adhesive cells have a tendency 
to increase their contact interfaces [22] , not to minimize them. 
Lecuit and Lenne [14] recently reviewed a large number of ex- 
periments, and show that a cell's surface tension results from 
the opposite actions of adhesion and cytoskeletal contraction. 
These are the ingredients of a second model [231 124] . 

Our approach is to find out if the observed cell packings 
and shapes can be described with one of these models, based 
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on the knowledge we have from the experiments. With min- 
imal and realistic assumptions, only the second model repro- 
duces the topology and geometry of the wildtype and mutant 
ommatidia. 

This shows that the competition between adhesion and 
cell cortex tension is needed to describe this specific cell pat- 
tern. We thus confirm and refine the conclusion that surface 
mechanics are involved in the establishment of cell topology 
and geometry. Adhesion plays an important role therein, but 
its role can only be understood when taking into account its 
effect on the cortical cytoskeleton. 

Results 

Model simulations. Why certain shapes are observed more of- 
ten than others depends on the developmental history of the 
tissue, which is determined by e.g. the sequence of cell dif- 
ferentiation, and cell divisions and deaths. Since a lot is still 
unknown about the developmental history, we do not include 
it in the modeling. However, since cells seem in mechanical 
equilibrium at any moment in development (c.f. 1141). future 
insights in developmental gene regulation could be translated 
in parameter changes that permit the modeling of the dynam- 
ics of development. 

Simulations thus start from unstable initial conditions (SI 
Fig. [8]l designed to favor the random search of final, stable 
topologies. We do not expect to find a quantitative corre- 
spondence between the frequency of topologies in simulations 
and experiments. We regard only the final result of the model 
simulations: we have found a local equilibrium, when the sim- 
ulated shape does not change anymore. 

We compare this shape with the experimental results 
(topology, geometry). Distinguishing between topologies is 
trivial. But, due to the variability of membrane fluctuations, 
we find that it is difficult to describe the geometrical char- 
acteristics (e.g. contact angles for the mutant ommatidia, 
interface lengths, elongation of cells) by quantitative mea- 
surements: one obtains more information by looking at the 
image ('eyeballing'). Quantitative measurements serve as a 
complement to the eyeballing when enough data are available 
(c.f. Fig. [TJ, not as replacement. We determine for each 
model which parameters do influence the shape of the cone 
cells; for the other parameters, we choose reasonable values 
(e.g., a compromise between simulation speed and precision, 

c.f. US]). 

We assume that (i) the adhesion strength is determined 
by the presence of these cadherins: when the two of them 
are present (i.e. at interfaces between cone cells), adhesion is 
thus stronger. Mutants should be modeled by only changing 
existing parameters. We thus require that (ii) to model the 
7?oi-mutants, we only need to change the number of cone cells; 
(iii) to model the cadherin mutant ommatidia, only the adhe- 
sion for the mutant cells should be changed (i.e. diminished 
for deletion, increased for overexpression) ; and (iv) all cells of 
a cell type that share the same mutation should be modeled 
using the same parameter values. 

Constant tension model. A stronger adhesion between cells i 
and j is represented by a lower interfacial tension 1 1 31 1171 1151 
116] . "/ij > 0, which is a constant depending only on the cell 
types of i and j. We minimise the energy: 

£= Yl T^j+Aa^A-Ak) 2 • [1] 

interfaces cells 



Pij is the length of the interface between cells i and j, Ai is 
the cell's area (the 2D equivalent of volume), Aoi is the cell's 
preferred area (target area), and Xa is the area modulus (a 
lower value allows more deviations from Ao). The values of 
Aoi are inferred from the experimental pictures, with cone (C) 
cells being smaller than primary pigment (P) cells. 

We assume C-C adhesion "foe, mediated by both E- and 
N-cadherin, to be stronger than C-P and P-P adhesion, yep 
and 7pp, which are mediated by E-cadherin alone. We as- 
sume the latter two to be equal: jee < Jgp — 7pp- Only 
three parameters need to be explored extensively: 7cc, 7c p 
(= 7pp), and Xa- The tensions 7 influence the cell shapes 
directly, whereas Xa determines a cell's deviations from the 
target area. 

Starting the simulations with a four-cell vertex (SI Fig. 
[8jA) , we systematically find an incorrect topology (Fig. [2j\): 
the anterior and posterior C cell touch. Even if we force the 
correct one, where the polar and equatorial C cell touch, it 
is unstable and decays into the incorrect one: the interfaces 
between the P cells are under tension, and pull the polar and 
equatorial C cells apart. 

To obtain the correct topology, we need another assump- 
tion: either that the adhesion between polar and equatorial 
C cells is stronger (Fig. [2J3) ; or that the P cells pull less on 
them (by having a stronger adhesion, Fig. [2p). Still, the ge- 
ometry is quite different from the experiments: notably the 
interface between the polar and equatorial C cell is too short 
in simulations. Besides, there is no experimental evidence to 
support these assumptions. 

Another optimization strategy is to determine (up to a 
prefactor) the tensions of three interfaces A, B, D that meet in 
a vertex from the experimentally observed contact angles a, f3, 
8 (a + /3 + 5 = 360°) by using 7a/ sin a = jp/ sin/3 = 70/ sin 5 
|26l 112] . We inject those tensions in the model. By construc- 
tion we obtain the correct contact angles, and thus topology; 
but the overall geometry (especially the interface lengths) dif- 
fers considerably from observations (results not shown). 

For the mutant ommatidia, the requirements (ii) to (iv) 
mentioned above could not be satisfied with this model: there 
are too many cases where other parameters need to be changed 
as well. We conclude that this model is insufficient to coher- 
ently describe the experiments. 

To obtain the observed shapes, it would certainly be pos- 
sible to choose a tension for each individual interface. But 
if the tension was just an input parameter without biological 
basis, then the model would not be predictive, nor help to 
understand the differences between the cells. 

Variable tension model. Adhesion between two cells tends to 
extend their contact length; it thus contributes negatively to 
the energy, —JijPij, where J > 0: in agreement with intu- 
ition, a higher J describes a stronger adhesion, while J — 
in absence of adhesion [231 124] . 

This extension is compensated by an elastic cell cortex 
term, Xp (Pi — Poi) 2 , where Xp is the perimeter modulus, and 
Poi is the target perimeter of cell i. The cell perimeter is the 
sum of its interfaces, Pi — J^Pij. We thus minimise the 
energy: 

£ = - J v p v+J2 [M^-PoO'+Aa^-AoO 2 ] • 

interfaces cells 

[2] 
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The interfacial tension 7^ = d£ij /dPij between cells i and 
j is the energy change associated with a change in membrane 
length (c.f. pi]); eq. ^ yields: 



Hi 



-Jy + 2A P (Pi - Poi) + 2A P (P, - P 



[3] 



As in the previous model, 7y is positive, else the cell would 
be unstable. However, it is no longer an input parameter. 
A stronger adhesion (high J) decreases the tension: this will 
usually cause an extension of the perimeter, which increases 
this tension even more. 

We represent all adhesion terms as combinations of E- and 
N-cadherin mediated adhesion (Je and Jn, respectively). In 
the wildtype, the adhesion between C cells is mediated by 
both cadherins, so Jcc = Je + Jn', whereas all other inter- 
faces only have E-cadherin, so Jpp = Jcp = Je- Values 
of Ao are estimated from pictures. The target perimeter Po 
(expressed in units of 2\/'kAo) should be larger for cells that 
deviate more from a circular shape, i.e. for the P cells. 

We thus adjust 6 main parameters: Je, Jn, Poc, Pop, 
Xp, Xa, which is too much to explore systematically. We 
adjust the parameters by hand, for wildtype and mutant con- 
figurations simultaneously, since the wildtype alone does not 
sufficiently constrain the number of optimal parameter com- 
binations. 

Unless indicated, throughout this paper, and for all figures 
except Fig. [2] we use eq. Q with the same set of parameters 
(SI Table [I]) for wildtype and mutants. 

Wildtype. Starting the simulations with a four-cell vertex (SI 
Fig. |8j\), the cells relax either into the correct topology where 
the polar and equatorial cells touch (Fig. Ilj; or into the incor- 
rect one where anterior and posterior cells touch (analogous 
to Fig. [5J\). Both topologies are stable, i.e. they are local 
energy minima. 

In the correct topology, the geometry of the simulated 
ommatidium resembles well the experimental pictures. More 
quantitatively, the contact angles measured in simulations and 
in experiments agree as well (Fig. [T]). In contrast to the con- 
stant tension model, we do not need additional assumptions. 

We found that the adhesion of secondary and tertiary pig- 
ment cells should be much stronger than can be expected from 
E-cadherin alone (J23 > Je, SI Table [TJ, otherwise they loose 
contact. Experimentally, deleting the E-cadherin of these 
cells does not induce any geometrical or topological change 
[9]. Both experiments and simulations thus suggest that sec- 
ondary and tertiary pigment cells might have other adhesion 
molecules than E- and N-cadherin. 

i?02-mutants. Without any additional parameter, we can sim- 
ulate different numbers of C cells (Poi-mutants); the total size 
of the simulation lattice is adjusted accordingly. For one, two, 
three and five C cells, only one topology is observed in exper- 
iments, and the same one in simulations (SI Fig. j9|. 

For six C cells, three topologies are observed experimen- 
tally (Fig. [3j\-C). Theoretically, there are two more possible 
equilibrium topologies for 6-cell aggregates, which are never 
observed although one of them has a smaller total interface 
length (simulations using the Surface Evolver, S. Cox, unpub- 
lished results 2004). We here performed a total of 42 Potts 
model simulations with different random seeds (see Methods), 
and found only three topologies (Fig. [3JD-F): they correspond 
to the observed ones. 



We observe in Fig. [3J\ and C that the entire ommatidium 
is elongated. Besides, ommatidia of Poi-mutants do not all 
have six sides and are assembled into a disordered pattern (see 
[5]). Thus, in Poi-mutants, ommatidia have variable shapes, 
which origin is not easily understood (especially for mutants 
with more pigment cells). Since in turn the shape of the om- 
matidium influences the geometry of its C cells (results not 
shown), studying the geometry of the C cells in more details 
would only by possible by adding more free parameters. 

N-cadherin mutants. Again without any additional parame- 
ter, simply by suppressing Jn, we could predict the pattern 
of ommatidia with N-cadherin deficient C cells. Since N- 
cadherin is only present on interfaces between C cells, deletion 
means we set the adhesion between mutant and wildtype C 
cells as J c c = Jcc = Je (mutant cells are denoted by lower 
case letters). 

We predict the correct topologies (Fig. [4|\-F and I-N), 
most of which are the same as in the wildtype. We predict 
qualitatively the main geometrical differences between mu- 
tants and wildtype: (i) the length of the interfaces between 
mutant cells and wildtype C cells decreases; (ii) the contact 
angles change; (iii) the interface length between the remaining 
wildtype C cells increases (Fig. |4|\-B and I- J); and (iv) the 
length of the central interface increases (Fig. |4p and L). 

When the polar or equatorial cell is the only C cell with- 
out N-cadherin, we simulate (Fig. [4jVI-N) both topologies that 
coexist in experiments (Fig. |4j3-F). 

To simulate one mutant P cell that mis-expresses N- 
cadherin, we optimize Jc P - While for the wildtype Jcp = 
Je = 150, we find an increase for the mutant, Jc P = 150+600. 
The high adhesion of this P cell with the C cells severely dis- 
rupts the normal configuration. Many topologies that dif- 
fer considerably form the wildtype are observed in exper- 
iments and simulations (e.g. Fig. [4ji,P). When both P 
cells mis-express N-cadherin, they balance each other and the 
topology is back to normal (Fig. [4j31). Optimization yields 
J pp = 150 + 700 (Fig. \4p and SI Fig. [ToJ). Both J Cp 
and Jpp are higher than the wildtype value of C-C adhesion 
[Jcc = Je + Jn = 150 + 450). 

E-cadherin mutants. The mutant C cell in Fig. [5}\ does not 
express E-cadherin, and it lacks adherens junctions at the in- 
terfaces with the P cells [9]. To simulate it, it would seem 
natural to suppress Je at all interfaces, that is, J c p = and 
JcC = Jn- With this assumption, we obtain the correct topol- 
ogy, which is the same as in the wildtype; however, the simu- 
lated geometry (not shown) is also the same as the wildtype, 
while the experiment is significantly different (Fig. [5|\). If we 
rather assume that C-C adhesion is unchanged by this muta- 
tion (JcC = Jcc), we obtain a good agreement (Fig. [SjD). 

E-cadherin overexpression in C cells (but not in P cells) 
significantly affects the pattern, yielding a coexistence of dif- 
ferent topologies: in Fig. [6j\ and B, the same cells are mu- 
tants, but the topologies differ; the same holds for Fig. [6] D 
and E. We predict the observed topologies (all stable) and, 
qualitatively, the geometries (Fig. [6|?-J) when we increase 
the C-P cell adhesion from Jcp = Je — 150 to J c p = 300; 
while we find that the adhesion between a wildtype and mu- 
tant C cells should not change, J c c = Jcc = Je + Jn, we 
should change it if both are mutants, J cc = 350 + Jn- Since 
E-cadherin overexpression in P cells rarely induces geometri- 
cal or topological changes [9] , we do not change their adhesion 



Footline Author 



PNAS I Issue Date | Volume | Issue Number | 3 



values. 

E- and N-cadherin mutants. We predict the effect of both 
E-cadherin and N-cadherin missing in C cells by setting 
J c c = JcP = Jcc = 0. Mutant C cells do not adhere to 
any of their neighbors, Fig. [5J3-F: intercellular space becomes 
visible between the cells, and the cells have shrunken. This 
agrees well with experiments, where mutant C cells lose the 
apical contacts with their neighbors (Fig. [HJ3-C). 

Discussion 

Constant and variable tension models. When surface tension 
is a constant model parameter, only modified by adhesion, the 
surface mechanics are soap-bubble-like: minimization of the 
interfaces with cell type dependent weights |13U17l[T51ll6H27| . 
This model proves to be insufficient here. However, in studies 
focussing on larger aggregates (10 2 to 10 4 cells) |17l 1271 [28] . 
constant surface tension was sufficient to explain tissue round- 
ing and cell sorting, and even Dictyostelium morphogenesis 
|29| . This constant tension model catches two important fea- 
tures of tissues of adherent cells: first, cells tile the space 
without gaps or overlap; second, the interface between cells is 
under (positive) tension, which implies for instance that three- 
cell vertices are stable, unlike four-cell ones, [101 111] and thus 
severely constrains the possible topologies |llj . 

In the present example of retina development, we show 
that interfacial tension should be variable, as described in a 
second model 23 , 24 . Tension results from a adhesion-driven 
extension of cell-cell interfaces, balanced by an even larger 
cortical tension (eq. |J3|). It explains correctly the topologies 
of many observations, and correctly simulates the geometries. 
It requires more free parameters; but they are tested against 
many more experimental data; and their origins, signs and 
variations are biologically relevant |14| . 

Adding more refinements (and thus more free parameters) 
would be easy, but does not seem necessary to describe the 
equilibrium shape of ommatidial C cells. The parameters 
should not be taken as quantitative predictions, since in vivo 
biophysical measurements to compare them to are lacking. 

Adhesion. By adjusting a set of 6 independent free parame- 
ters in this variable tension model, we obtain topological and 
geometrical agreement between the simulations and the pic- 
tures of 16 different situations: the wildtype (Fig. [TJ, the six 
topologies observed in the 7?oi-mutants (Fig. [3] and SI Fig. 
9 1; as well as the nine cadherin deletion mutants (Figs. [4]A_-F, 
5 1 by setting the corresponding parameter to zero. 

We also simulate 7 cadherin overexpression mutants, by 
re-adjusting the corresponding parameter (Figs. [4]G!,H, [6|: 
adhesion is increased. The strongest increases are found when 
two overexpressing cells touch: this corresponds to the idea 
that the adhesion strength depends on the availability of cad- 
herin molecules in both adhering cells. 

We found two cases where a mutation does not seem to 
change the adhesion strength: first, when deleting E-cadherin 
from one C cell, its adhesion with a normal C cell is un- 
changed (Fig. [5j3); second, we rarely observed shape changes 
in E-cadherin overexpressing P cells in experiments (c.f. [9]). 

Indeed, while a linear relation between cadherin expres- 
sion and adhesion strength has been found in vitro [30] . this 
need not be true in vivo, since cells have many more ways to 
regulate protein levels. These exceptions, thus, do not con- 



tradict the conclusion that the shapes observed in mutants 
are the effect of altered adhesion: an increase in the case of 
overexpression, a decrease in the case of deletion. 
Cortical tension. In the variable tension model, the perime- 
ter modulus \p and the target perimeter Po reflect the role 
of the cortical cytoskeleton. The target perimeter is always 
smaller than the perimeter, therefore the interfacial tension 
7ij (eq. [J3j| ) is always positive, else the cell would be unstable 
and fall apart or disappear. The cortex of the simulated cells 
is contractile, and generates tension. This tension depends 
on the perimeter P of the cell, which length depends on the 
cell's shape, which in turn depends on the tension: there is a 
feedback between tension and shape, and thus between each 
cell and its neighbours. 

To understand the effect of this feedback, let us consider 
the wildtype ommatidium. We assume that the four C cells 
have equal adhesion properties. The tension at the interfaces 
between the two P cells pulls at the polar and equatorial C 
cell. When the tension is constant, these cells will therefore 
be pulled apart (Fig. |2|\.): the cells do not react on their de- 
formation. When the tension, however, depends on the cell's 
perimeter, pulling at those cells deforms them, and increases 
their tension: energy minimization thus requires that they 
stay in contact. 

The prediction that cytoskeletal contractility is essential 
for the establishment of cell shape should be tested, e.g. by 
treating the cells with cytoskeletal inhibitors |19l 131] . or ge- 
netically modifying the cytoskeleton. Since the cytoskeleton 
has multiple functions that could interfere with adhesion (c.f. 
[321 15]). the results will be difficult to interpret. Preliminary 
experimental results (not shown) do indicate that genetically 
disturbing Rho-family GTPases influences the cell shape. The 
role of the cytoskeleton has been confirmed in various tissues 
and organisms (see |14l 133] for reviews). We here present a 
computational framework able to test this hypothesis, which 
can be extended to other tissues, ranging from patterns of few 
cells to large-scale aggregates. 

Methods 

Experiments. Retinas were stained and analyzed as described 
in [5] and [34] . In short, cell contours were visualized by stain- 
ing either with cobalt sulfide (Fig. [3] and SI Fig. [9b, or with 
the antibodies against DE-cadherin, DN-cadherin (referred to 
as E- resp. N-cadherin in the rest of the text), /3-catenin 
or /3-spectrin. Rough eye (Roi) flies were used to examine 
the topology and geometry of variable number of cone cells. 
The effect of eliminating or overexpressing cadherin molecules 
was studied in mosaic retinas composed of wild-type and mu- 
tant cells (see [9])- We examined more than five retinas in 
each experiment. Thus at least several hundreds of omma- 
tidia (> 500) were examined for the wildtype and each muta- 
tion, except the E- and N-cadherin overexpression, in which 
case approximately 100 ommatidia were examined each. Some 
pictures used for the analysis were published previously [9ll34|. 
Model simulations. The cellular Potts model [l5l[T6| is a stan- 
dard algorithm to simulate variable cell shape, size and pack- 
ing [25]. Its use in biology is motivated by the capability 
to handle irregular, fluctuating interfaces (c.f. [35]); the pix- 
elisation induced by the calculation lattice can be chosen to 
correspond to the pixelisation in the experimental images. 
Each cell is defined as a certain set of pixels, here on a 2D 
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square lattice; their number defines the cell area A. The cell 
shapes change when one pixel is attributed to one cell instead 
of another. Our field of simulation for one ommatidium is a 
hexagon with sides of approximately fOO pixels (its surface 
is Ahex = 25160 pixels, about the same as in experimental 
pictures). We use periodic boundary conditions, as if we were 
simulating an infinite retina with identical ommatidia. Ini- 
tially, the whole hexagon is filled with cells, approximately at 
the right positions (SI Fig. |s|) - We treat bristle cells as tertiary 
pigment cells: both are situated at the edge of three omma- 
tidia. These initial conditions, with an unstable n-cell vertex 
in the middle, do not fix the final configuration in advance. 
Simulations can be started with different seeds of the random 
number generator, to explore whether multiple solutions are 
possible. 

Shape is relaxed in order to decrease the energy £ , eq. ijlj 
|15| or eq. Q [53]. The algorithm to minimize £ uses Monte 
Carlo sampling and the Metropolis algorithm, as follows. We 
randomly draw (without replacement) a lattice pixel, and one 
of its eight neighboring pixels. If both pixels belong different 
cells, we try to copy the state of the neighboring pixel to the 
first one. If the copying diminishes £ , we accept it; and if it 
increases £ , we accept it with probability P — exp (— A£/T). 
Here A£ is the difference in £ before and after the considered 
copying. The prefactor T is a fluctuation (random copying) 
allowance: it determines the extent of energy-increasing copy 
events, leading to membrane fluctuations 35:. Since all energy 
parameters are scalable with the fluctuation allowance T, we 
can fix it without loss of generality; for numerical convenience 
we choose numbers of order of a hundred. 

We define one Monte Carlo time step (MCS) as the num- 
ber of random drawings equal to the number of lattice pixels. 
It takes approximately 600 to 4000 MCS to attain a shape 
that does not evolve anymore, that is, in mechanical equi- 
librium where stresses are balanced. We run the simulation 
much longer (up to 10 6 MCS) to test if topological changes 
occur. 

To avoid possible effects of lattice anisotropy on cell 
shapes, we compute P and £ by including interactions up 
to the 20 next next nearest neighbours [36]. All perimeters 
indicated here are corrected by a suitable prefactor 10.6 to 
ensure that a circle with an area of A pixels has a perimeter 

2V^a [37]. 

In experiments, interstitial fluid is present in small 



amount, and cells can lose contact (Fig. J5]B,C). To simu- 
late it in our 2D model, at each MCS we randomly choose one 
pixel at a cell interface and change its state into 'intercellu- 
lar space' (a state without adhesion, nor area and perimeter 
constraints). In addition, we choose the sum of all cells target 
areas to be less than the total size of the hexagonal simu- 
lation field (Eceiis a <k = 0-95 A hex , SI Table [B. Only when 
cells lose adhesion (J = 0) do we actually observe intercellular 
space in simulations (Fig. [5}3,F). 

We try different parameters and adjust them to improve 
the visual agreement ('eyeballing') between simulated and ex- 
perimental pictures. To estimate our uncertainty, we note 
that 5 — 10% changes in the values of the adhesion parameters 
do not yield visible changes in the geometry, while 10 — 30% 
changes do; see SI Fig. [l0]for an example of the determination 
of Jc P - 

Images. Once we simulate the correct topology, we measure 
the contact angles of straight lines fitted through the inter- 
faces that meet in the vertex. The line should be long enough 
to avoid grid effects; we fit a straight line using the first 15 
first-order neighboring sites. Since the simulated cells show 
random fluctuations, statistics are obtained by measuring the 
contact angles several times during the simulation, or in sim- 
ulations with different random number seeds. 

In experimental pictures, we measure contact angles in 22 
wildtype ommatidia by hand, aided by the program ImageJ 
|38] . Ommatidia have two axes of symmetry, and we consider 
the ommatidia to consist of four equal quarters, which gives 
us 88 measurements for each angle (and 44 measurements of 
the angles that are intersected by the axes of symmetry). The 
variation between different wildtype ommatidia is larger than 
in simulations (Fig. [T]). In mutant ommatidia, the error bar 
is even larger, so that we did not attempt at any quantitative 
comparison. 
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measured angle in simulation (degrees) 



Fig. 1: Wildtype. Contact angles measured in experiments and in simulations are plotted as average ± statistical standard 
deviation; the straight line represents y — x. Inset left: an ommatidium stained for E-cadherin; anterior (a), posterior (p), 
polar (pi) and equatorial (e) cone cells. Inset right: variable tension model simulation, with cone cells (C), and primary (P), 
secondary (2) and tertiary (3) pigment cells. One ommatidium contains four times the angles a, /3, 5, rj, and (, and two times 
e and 8. 

ww<s> 

Fig. 2: Constant tension model simulations. A: 'ycc ~ 40, 7cp = 7pp = 80. B: Same as (A), but with lower tension (stronger 
adhesion) between the polar and equatorial cone cell, 7 po i ar , equatorial = 20. C: Same as (A), but with lower tension (stronger 
adhesion) between the primary pigment cells: 7pp = 40. 




Fig. 3: Roi mutants with 6 cone cells. A-C: Experimental pictures, from ref. |S]. D-F: Corresponding simulations. 




Fig. 4: N-cadherin mutants. Mutant cells are indicated with a "+" for overexpression, "-" for deletion. A-H: Experimental 
pictures, A-D from ref. [9]. I-P: Corresponding simulations. 
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Fig. 5: Loss of adhesion. Mutant cells are indicated with a A: A mutant cone cell lacking E-cadherin. B-C: Double 
mutant cone cells for E-cadherin and N-cadherin. D-F: Corresponding simulations. A-C from ref. [9]. 
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SI Fig. 7: Confocal microscopy images of the Drosophila retina. The pictures measure 100 /im x 100 /im. a) /3-catenin, a 
component of the adherens junction, is stained green. Nearly all catenin fluorescence between the cone cells is seen in a layer of 
1.26 fim thick. One ommatidium consists of four cone cells (c) and two primary pigment cells (p), surrounded by six secondary 
(2) and three tertiary (3) pigment cells and three bristle cells (b). In this particular ommatidium, one bristle cell is replaced by 
a tertiary pigment cell. The cone cells can be subdivided into a polar (pi), equatorial (eq), anterior (a) and posterior (po) cone 
cell, according to their position, b) N-cadherin fluorescence in the same plane of focus. N-cadherin is restricted to the cone - 
cone interfaces. 



SI Fig. 8: Initial conditions of each simulation with four cone cells and two primary pigment cells (left), and six cone cells and 
three primary pigment cells (right). Periodic boundary conditions imply that the secondary pigment cells (purple) and tertiary 
pigment cells (red) that are marked with the same symbol, are treated as parts of the same cell. 



SI Fig. 9: Experiments and simulations showing ommatidia with 2 (A,D), 3 (B,E) and 5 (C,F) cone cells. A-C from ref. [9]. 



SI Fig. 10: Determination of the adhesion between cone cells and two N-cadherin mis-expressing pigment cells. Simulations 
are shown with values J Cp = 150 (A), 600 (B), 700 (C), 750 (D), 800 (E), 850 (F). (A) corresponds to wildtype, (D) corresponds 
best to the mis-expression experiment (Fig. Hp). 



SI Table 1: Simulation parameters of the wildtype in the variable tension model, Free parameter adjusted to compare to 
wildtype observation, t: Parameters which value has little effect on the images. J: Parameter of no effect on the images, 
since cone cells almost never touch secondary or tertiary pigment cells. Target perimeters are expressed as a factor times the 
perimeter of a circle having the specific target area. E.g. a prefactor of 1 indicates that the target perimeter of the cell equals 
the perimeter if the cell is round and has an area equaling its target area. A cell with a prefactor > 1 (like the primary pigment 
cells) can deviate much from a round shape. Abbreviations: N, N-cadherin; E, E-cadherin; C, cone cell; P, primary pigment 
cell; 2, secondary pigment cell; 3, tertiary pigment cell. 
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